Typeset with jpsj2.cls <ver.l.2> 



Full Paper 



O 
O 



Polarized Neutron Inelastic Scattering Study of the Anisotropic Magnetic 
Fluctuations in the Quasi-ID Ising-like Antiferromagnet TIC0CI3 

Akira OOSAWA*, Kazuhisa Kakurai 1 , Yoichi Nishiwaki 2 and Tetsuya Kato 3 

Department of Physics, Sophia University, 7-1 Kioi-cho, Chiyoda-ku, Tokyo 102-8554, Japan 
1 Quantum Beam Science Directorate, Japan Atomic Energy Agency, Tokai, Ibaraki 319-1195, Japan 



'Department of Physics, Tokyo Institute of Technology, Oh-okayama, Meguro-ku, Tokyo 152-8551, Japan 
3 Faculty of Education, Chiba University, 1-33 Yayoi-cho, Inage-ku, Chiba 273-8522, Japan 

(Received February 3, 2008) 

Polarized neutron inelastic scattering experiments have been carried out in the quasi- 
1D Ising-like antiferromagnet TIC0CI3. We observed the longitudinal magnetic fluctuation 
S zz (Q,ui) for the spin- wave excitation continuum, which has not been observed in the un- 
polarized neutron inelastic scattering experiments of the quasi- ID Ising-like antiferromagnets 
CsCoCbj and TIC0CI3 so far, together with the transverse magnetic fluctuation S xx (Q,uo). We 
compared both obtained intensities of S xx {Q,oj) and S zz (Q,ui) with the perturbation theory 
from the pure Ising limit by Ishimura and Shiba (IS), and a semi-quantitative agreement was 
found. 
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1. Introduction 

The magnetic excitations of the ID Ising-like antifer- 
romagnet with the following Hamiltonian, 



H = 2J Y1 i S j S Ui 



( QX QX 



(1) 



have been attracting much attention because the exci- 
tations are characterized as the domain-wall excitations 
and the behavior of the excitations is much different 
from that of the conventional spin-wave excitations of 
the Heisenberg antiferromagnet. 

Previously, the unpolarized neutron inelastic scatter- 
ing experiments have been performed on the quasi-lD 
Ising-like antiferromagnet CSC0CI3. 1,2 In the experi- 
ments, two magnetic excitations, which arc the cen- 
tral mode corresponding to the motion of thermally ac- 
tivated domain walls firstly predicted by Villain 3 and 
the spin-wave excitation continuum corresoponding to 
the domain-wall pair excitation, were observed. For the 
central mode, both the transverse magnetic fluctua- 
tion S xx (Q,ui) and the longitudinal magnetic fluctua- 
tion S zz (Q,uj) were observed and the separation of the 
two fluctuations has been successfully performed using 
the property of the neutron magnetic inelastic scattering 
cross section, namely the magnetic fluctuations perpen- 
dicular to the scattering vector Q is only contributing 
to the cross section. While, for the spin-wave excitation 
continuum, the separation of both fluctuations failed be- 
cause the intensity of S zz (Q, to) is much smaller than that 
of S xx (Q,u>) so that S zz (Q,uj) is masked by S xx (Q,uj) 
out of the scattering plane which is always observed in 
the unpolarized experiments. 

In this paper, we demonstrate the advantages of the 
polarized neutron inelastic scattering experiments for the 
investigation of the anisotropic magnetic fluctuations in 
the ID Ising-like antiferromagnet, which enables the sep- 



aration of S xx {Q,uj) and S zz (Q,ui), and report the po- 
larized neutron inelastic scattering investigation on the 
quasi-lD Ising-like antiferromagnet TIC0CI3. 

TIC0CI3 has the CsNiCi3-type hexagonal crystal struc- 
ture (space group symmetry P&^/mmc) at room tem- 
perature, in which the magnetic Co ions form the linear 
chain along the c-axis and these chains make triangular 
lattice in the c-plane. Because Co ions have a large Ising 
anisotropy along the c-axis due to the Kramers doublet, 
this system can be expressed as the quasi- ID Ising-like 
antiferromagnet. From the dielectric constant measure- 
ments, it was found that this system undergoes the suc- 
cessive structural phase transitions with fcrroelectricity 
at T st2 = 165 K, T st3 = 75 K and T st4 = 68 K. 4 Besides, 
it was found from the magnetization measurements that 
the magnetic phase transition occurs at Tn=29.5 K. 4 The 
neutron diffraction measurements have been carried out 
in this system in order to investigate the crystal and 
magnetic structures. It was found that the crystal struc- 
ture varies from the hexagonal structure {PQ^/mmc) at 
room tempearture to the orthorombic one (Pbca) be- 
low T st 4 and the magnetic Bragg reflections indicative 
of the up-up-down-down-type magnetic structure have 
been observed below Tn at Q = ^, I) with odd h and 
odd I, which are assigned by the indices of the room- 
temperature hexagonal lattice. 5 The unpolarized neu- 
tron inelastic scattering measurements have been per- 
formed in order to investigate the magnetic excitations in 
TIC0CI3. 6 The spin- wave excitation continuum was ob- 
served. According to the IS theory, 7 we analyzed the ob- 
served spin-wave excitation, and the exchange constant 
2 J and the anistropy e were estimated as 14.7 meV and 
0.14 in TIC0CI3, respectively 

2. Experimental Details 

Polarized neutron inelastic scattering experiments us- 
ing the uniaxial polarized neutron scattering technique 8 
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Fig. 1. Scattering plane of the present experiment. The gray ar- 
rows denote the scattering vectors Q in the present measure- 
ments. 8 is the angle between the scattering vector Q and 
(h, h, 0)-axis. The magnetic fluctuations along each directions 
arc also shown. 



were performed using JAEA-TAS1 spectrometer in- 
stalled at JRR-3M in Tokai. The constant-fcf mode 
was taken with a fixed final neutron energy Ef of 
14.7 meV. Heusler(lll)-Hcusler(lll) monochromator- 
analyser configuration was used in the present experi- 
ments. The polarization was measured on the nuclear 
Bragg reflection and was determined to be around 96 %. 
A pyrolitic graphite-filter was used to suppress the higher 
order contaminations. The collimations were set as open- 
80'-open-open in order to gain intensity. We used a sin- 
gle crystal of TIC0CI3 with a volume of approximately 
0.5 cm 3 . The present measurements were performed at 
T = 34 K, which is just above Tn- The crystallo- 
graphic parameters were determined as a = 6.84 A and 
c = 5.93 A at T = 34 K, which arc defined in the room- 
temperature hexagonal lattice. The sample was mounted 
in the cryostat with its (h, h, 0)- and (0, 0, ^)-axes in 
the scattering plane, as shown in Fig. 1. The indices of 
the room-temperature hexagonal lattice are used in the 
present experiments. The relation between the indices 
of the room-temperature hexagonal lattice and those of 
the low-temperature orthogonal lattice is shown in ref. 
5. In the present experiments, the guide field was ap- 
plied perpendicular to the scattering (h,h,l) plane, i.e. 
vertical field configuration. In this configuration, the in- 
elastic magnetic scattering intensity in the spin-flip (SF) 
channel is proportional to the intensity of the magnetic 
fluctuations perpendicular to the vertical guide field, i.e. 
fluctuations in the scattering plane, while that of the 
non-spin-flip (NSF) channel is proportional to the inten- 
sity of the magnetic fluctuations parallel to the guide 
field, namely fluctuations out of the scattering plane as 
follows: 



-^NSF 

Isf 



C-S xx {Q,u>) 



(2) 



= C (sin 2 9 ■ S XX {Q, u) + cos 2 9 ■ S zz (Q, w)) 
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Fig. 2. Profiles of the constant-Q energy scans in TIC0CI3 for 
NSF and SF scatterings at Q = (0, 0, 2.5). The solid and dashed 
lines denote the fits by a Gaussian function and the energy de- 
pendent backgrounds, respectively. The horizontal bar indicates 
the calculated resolution width at E = 14.7 meV. 



where C and 9 are the scale factor and the angle be- 
tween the scattering vector Q and (h, h, 0)-axis, respec- 
tively. Hence the longitudinal fluctuation S zz (Q,u>) can 
be separated from the transverse fluctuation S xx (Q,lu) 
by measurements at a single (Q, uS) point. 

3. Results and Discussion 

First, we measured the NSF and SF scattering inten- 
sities for Q = (0,0,2.5) (Point-A in Fig. 1 ). Because 9 
is equal to 90° at Q = (0, 0, 2.5), we can expect that the 
same intensities, i.e. C ■ S xx (Q,u>), are observed for both 
NSF and SF scatterings. Figure 2 shows the profiles of 
the constant-Q energy scans in TIC0CI3 for NSF and SF 
scatterings at Q = (0, 0, 2.5). As shown in Fig. 2, the 
excitations with resolution limited and same intensities 
are clearly seen for both scatterings. From the IS theory, 
it was found that S xx (Q,to) in the ID Ising-likc antifcr- 
romagnet can be expressed as the following expression 
for the peak position 



2 J(l - 8e 2 cos 2 tt/) 



(4) 



and the boundaries of the continuum extending from uj 



to cj+, given by 



2J(1 ±2ecos7rZ). 



(5) 



(3) 



Hence, at Q = (0, 0, 2.5), the continuum should be closed 
so that the magnetic excitation with resolution limited 
should be observed at lu = 2 J. Because the exchange con- 
stant 2 J has been already estimated as 14.7 meV in this 
system, 6 the observed excitations at Q = (0, 0, 2.5) are 
consistent with the IS theory. Note that this system un- 
dergoes the successive structural phase transitions from 
the hexagonal crystal structure to orthorombic one. 4 ' 5 
However, it was found that structural phase transitions 
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Fig. 3. Profiles of the constant-Q energy scans in TIC0CI3 for 
NSF and SF scatterings at Q = (1.5, 1.5, 0). The solid and dashed 
lines denote the guide for the eyes and the energy dependent 
backgrounds estimated from the scan profiles at Q = (0, 0, 2.5), 
as shown in Fig. 2, respectively. The horizontal bar indicates the 
calculated resolution width at E = 13.5 meV. 

scarcely affect S xx (Q,lu) because the magnetic excita- 
tions with same intensities were observed for both scat- 
terings. The sloping background in Fig. 2 can be ascribed 
to the transfer energy dependent counting time due to 
the chosen constant-fcf mode, where the incident neu- 
tron energy is varied while the outgoing neutron energy 
is fixed at Ef — 14.7 meV. Since we measure normalized 
on the incident beam monitor counts, the counting time 
increases as the transfer energy is increased and concomi- 
tantly the instrumental background increases. These Q- 
independent sloping backgrounds for both channels will 
be fixed in all the following profiles to compare the ob- 
served intensities with the calculation. 

Next, we measured the NSF and SF scattering inten- 
sities for Q = (1.5,1.5,0) (Point-B in Fig. 1 ). Because 
9 is equal to 0° at Q = (1.5, 1.5, 0), we can expect from 
eqs. (2) and (3) that only S xx (Q,uj) and S zz (Q,u) are 
observed for NSF and SF scatterings, respectively. Fig- 
ure 3 shows the profiles of the constant-Q energy scans 
in TIC0CI3 for NSF and SF scatterings at Q = (1.5, 1.5, 
0). As shown in Fig. 3, no excitation can be seen for 
the SF scattering while a little peak structure is seen 
for the NSF scattering. In the IS theory, the intensity of 
3 zz {Q,oS) can be expressed as 

where 

V = eJ(l + e- 2ml ) (7) 



Fig. 4. Profiles of the constant-Q energy scans in TIC0CI3 for 
NSF and SF scatterings at Q = (1.5, 1.5, 1). The solid and dashed 
lines denote the guides for the eyes and the energy dependent 
backgrounds estimated from the scan profiles at Q = (0, 0, 2.5), 
as shown in Fig. 2, respectively. The dash-dotted line denotes the 
NSF scattering intensity calculated from the Sxx (Q, u>) spectrum 
of the IS theory (eq. (9) ) with the appropriate parameters given 
in the text. The horizontal bar indicates the calculated resolution 
width at E = 13.5 meV. 

while that of S xx (Q,u>) can be expressed as 

[ ^^IVT - (w - 2J) 2 [1 - 2e cos ttZ 
S xx (Q,u)~< -^r 1 ] for \u-2J\<2\V\ 
[ otherwise. 

(8) 

As shown in eqs. (6) and (8), S zz (Q,u>) should become 
while S xx (Q,uj) is finite when \u> — 2J| < 2\V\ = 4eJ at 
Q = (1.5, 1.5, 0). This prediction is consistent with the 
experimental results, as shown in Fig. 3. 

In order to detect S zz (Q,oj), we measured the NSF and 
SF scattering intensities for Q = (1.5, 1.5, 1) (Point-C in 
Fig. 1 ). Because 9 is estimated as 21° at Q = (1.5, 1.5, 1), 
we can expect from eqs. (2) and (3) that the inten- 
sities corresponding to S xx (Q,u>) and 0.135 X2; (Q, u>) + 
0.87S zz (Q,uj) are observed for NSF and SF scatterings, 
respectively. Figure 4 shows the profiles of the constant- 
Q energy scans in TIC0CI3 for NSF and SF scatterings 
at Q = (1.5, 1.5, 1). In the NSF scattering, we observed 
the magnetic excitation with asymmetric spectral shape 
which is the same as that observed in the unpolarized 
experiments of CSC0CI3 1 ' 2 and TIC0CI3, 6 namely the 
excitation is attributed to the large transverse magnetic 
fluctuation S xx (Q,u>), while the spectral shape of the 
magnetic excitation observed in the SF scattering is dif- 
ferent from that in the NSF scattering. We can expect 
that this difference is due to the longitudinal magnetic 
fluctuation S ZZ (Q, u>). Note that this difference is not due 
to the phonon excitations because phonon excitations 
should be only observed in the NSF channel. 8 

In the following, we compare the obtained results at 
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Fig. 5. Energy dependence of the observed residual SF scatter- 
ing intensity corresponding to the pure S zz (Q,lu) contribution, 
as shown in eq. (13), at Q = (1.5, 1.5, 1). Here the energy depen- 
dent background is subtracted. In the estimation, we averaged 
the two neighboring data in order to reduce the statistical error. 
The dash-dotted line denotes the calculated residual SF scatter- 
ing intensity based on the S zz (Q,lj) spectrum of the IS theory 
without adjustable parameters (see text). The solid line is only 
a guide to the eyes. The horizontal bar indicates the calculated 
resolution width at E = 13.5 meV. 



Q = (1.5, 1-5,1) with the IS theory. First, we consider the 
NSF scattering intensity. As mentioned above, S XX (Q, uj) 
can be expressed as eq. (8). Hence, at Q = (1.5, 1.5, 1), 
the NSF scattering intensity can be expressed as follows: 



with 



S xx (Q,u) 



'NSF 



327re 2 J 2 



C • S xx (Q,lu) 



(9) 



16e 2 J 2 -(uj- 2J) [1 + 2e 
H] for \u-2J\<2\V\ 



otherwise 



(10) 

The NSF scattering intensity calculated using this 
S xx (Q,lu) with the exchange constant 2 J and the 
anistropy e of 14.7 meV and 0.14, respectively, is plot- 
ted in Fig. 4 by the dash-dotted line. In order to fit the 
NSF scattering intensity obtained in the present mea- 
surements, the scale factor C is estimated as 3000. As 
shown in Fig. 4, the spectral shape is not very well re- 
produced. Note that the convolution of the resolution 
function is not considered in this comparison. However, 
this difference is not due to the convolution because the 
observed spectrum is sharper than the calculated one. 
The SF scattering intensity at Q = (1.5, 1.5, 1) can be 
expressed as follows: 



Isf = C(0.13S xx (Q,lj) + 0.87S zz (Q,u>)) 



(11) 



uj - 2 J 
4eJ 



(12) 



and S xx (Q,uj) shown in eq. (10). Hence the residual SF 
scattering intensity corresponding to the pure S zz (Q,uj) 
contribution can be estimated from NSF and SF scatter- 
ing intensities, namely 



ires = Isf — 0.13/nsf 
= C -0.87S zz {Q,uj). 



(13) 



Note that there are no adjustable parameters in this 
comparison because the scale factor C is already esti- 
mated from the obtained NSF scattering intensity, as 
mentioned above. The observed and calculated residual 
SF scattering intensities after the subtraction of the lin- 
ear energy dependent background are shown in Fig. 5. 
For the estimation of I res , we averaged the two neigh- 
boring data in order to reduce the statistical error. This 
procedure can be considered to be valid because the ex- 
perimental resolution in the energy scan is approximately 
2.5 meV, as shown in Fig. 5. The intensity of the ob- 
served longitudinal fluctuation S zz (Q,uj) in the vicinity 
of E = 2 J = 14.7 meV is quite well reproduced by the IS 
theory. Note that the intensity of the longitudinal fluc- 
tuation S zz (Q,uj) is much smaller than the transverse 
fluctuation S xx (Q,uj) in accord with the theoretical pre- 
diction. The spectral shapes of the experimentally ob- 
served S zz (Q,uj) and S xx (Q,uj) do not quite well cor- 
respond to those calculated by the IS theory. There are 
some discussions, in which this difference is attributed to 
the additional terms, such as the exchange mixing 2 ' and 
the next-nearest-neighbor ferromagnetic interaction, 10, 11 
which are not considered in the IS theory. The compar- 
ison between the calculated spectra including the addi- 
tional terms and the observed ones in the present exper- 
iments is needed. 

4. Conclusion 

We have presented the results of the polarized neutron 
inelastic scattering measurements in the quasi-lD Ising- 
like antiferromagnet TIC0CI3. We clearly identified the 
longitudinal magnetic fluctuation S zz (Q,ui) for the spin- 
wave excitation continuum at w ~ 2 J, which has not 
been observed in the unpolarized neutron inelastic scat- 
tering experiments of the quasi-lD Ising-like antiferro- 
magnets CSC0CI3 and TIC0CI3 so far, together with the 
transverse magnetic fluctuation S xx (Q,uj). It was found 
that the observed longitudinal and transverse magnetic 
fluctuation spectra at characteristic wave vectors semi- 
quantitatively correspond to those calculated by the per- 
turbation theory from the pure Ising limit, as shown in 
Figs. 2 through 5. 
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